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We report that nonmagnetic heavy-fermion (HF) iron oxypnictide CeFePO with two-dimensional 
XY-type anisotropy shows a metamagnetic behavior at the metamagnetic field Hm — 4 T perpen- 
dicular to the c-axis and that a critical behavior is observed around Hyi. Although the magnetic 
character is entirely difi'erent from that in other Ce-based HF metamagnets, Hm in these metam- 
agnets is linearly proportional to the inverse of the effective mass, or to the temperature where the 
susceptibility shows a peak. This finding suggests that Hm is a magnetic field breaking the local 
Kondo singlet, and the critical behavior around Hm is driven by the Kondo breakdown accompanied 
by the Fermi-surface instability. 

PACS numbers: 76.60.-k, 71.27.-|-a 74.70.Xa 



Metamagnetism is represented by a sudden increase 
in magnetization with increasing an applied field. In 
heavy fermion (HF) systems, CeRu2Si2 with the tetrag- 
onal ThCr2Si2 structure shows the metamagnetic behav- 
ior at about 7.7 T when a magnetic field [H) is applied 
parallel to the c-axis. Although various experiments as 
well as theoretical studies have been carried outP,|3|, the 
mechanism is still controversial. In order to understand 
the metamagnetic behavior in HF systems, it might be 
desired to investigate new metamagnetic compounds. 

The iron oxypnictide CeFePO is a related material of 
the iron-based superconductor LaFeP0[l,|3]. They pos- 
sess the same two-dimensional layered structure, stacking 
Ce(La)0 and FeP layers alternately. Briining et al. re- 
ported that CeFePO is a magnetically nonordered HF 
metal with a Sommerfeld coefficient 7 = 700 mJ/(mol 
K^)0. At present, it is difficult to synthesize large 
single crystals of CeFePO for NMR measurements, but 
■^^P-NMR can probe in-plane and out-of-plane magnetic 
response separately using c-axis aligned polycrystalline 
samples. Here we report novel metamagnetic behavior 
observed in _L c, and suggest that metamagnetism of 
Ce-based HF compounds is driven by Kondo-breakdown 
(drastic reduction of c-f hybridization) as clarified ex- 
perimentally. 

The polycrystalline CeFePO was synthesized by solid- 
state reaction Basic properties are consistent with 
the previous report To measure anisotropic mag- 
netic properties of CeFePO, the samples were uniaxially 
aligned using a magnetic field 0. The polycrystalline Ce- 
FePO was ground into powder, mixed with stycast 1266, 
and was rotated in the external field of 1.4 T while the 
stycast cures. The c-axis of the sample is nicely aligned, 
which is shown from the angle dependence of ■^^P-NMR 
spectra (see in the inset of Fig[T]), and '^^P-NMR mea- 
surement was performed on the sample. 
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FIG. 1. (Color online) (Main panel) T dependence of //-swept 
NMR spectra at 31.7 MHz for 7/ _L c (solid line) and H \\ c 
(broken line). K = was determined by reference material 
H3PO4. (Inset) Angle dependence of //-swept NMR spectra at 
31.7 MHz measured at 20 K. is the angle between magnetic 
field and c-axis. Solid line is corresponding to fitting line. 

Figure [T] shows //-swept NMR spectra in // || c and 
// _L c obtained at 31.4 MHz and various temperatures 
(T). The resonance peak for // || c is almost T inde- 
pendent, but the peak for H J- c shows the character- 
istic T dependence originating from x(T). The Knight 
shift A'_L(||) was determined from the peak field of the "^^P 
NMR spectrum obtained in H perpendicular (parallel) to 
the c-axis. K ~ was determined by reference material 
H3PO4. Ki{T,H){i =_L and ||), which is the measure of 
the local susceptibility at the nuclear site, is defined as, 

i^.(T,//.e.)=f^^i-^) oC^^M^l^ (1) 

V /ires / ^ros 

where //res arc magnetic fields at resonance peaks, Hq 
and cjo are the resonance field and frequency of bare 
•^^P nucleus and have the relation of loq = jnHo with 
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FIG. 2. (Color online) (a) T dependence of the Knight shift 
determined at the peaks of 7? ± c and H \\C spectra obtained 
at various H. Strong anisotropy of Knight shift suggests that 
static spin properties possess XF-type spin anisotropy. (b) R 
dependence of magnetization Mi{H) {i =_L and ||) using the 
relation of Mi{H) = Ki{H)H^cs/ A]^{. Solid and broken lines 
are guide to eyes. M±{H) suddenly increases with increasing 
H and deviates from linear relation in the field range of 3 - 
5 T, which is a definition of a metamagnetic behavior, while 
such a behavior was not observed in My (H) up to 6.2 T and 
down to 1.5 K. 



gyromagnetic ratio 7„, and Mi{T,H^cs) is the magne- 
tization under i/i,ros(* =-L and ||) at T. A"|| is almost 
independent of T and iJ, whereas K± shows strong T 
dependence originating from the Curie- Weiss behavior of 
x(T) above 10 K as shown in Figllja). The anisotropic 
Knight shift suggests that static spin properties possess 
Xy-type spin anisotropy. It should be noted that K± 
exhibits H dependence below 4 K and above 2 T, indica- 
tive of a non- linear relation between and H . Us- 
ing the hyperfine coupling constant '^"'^^hf = 0.2 T/^b, 
which is estimated from the plot between isotropic com- 
ponent of K and x{T) above 10 K (not shown), we can 
plot Mi{H) against H in FiglJJb). M^{H) becomes su- 
perlinear against H at 0.1 K, which is the hallmarks of 
metamagnetism, whereas {H) is linear up to 6.2 T, 
which is again highly anisotropic. 

Next, we focus on T and H dependence of low-energy 
spin dynamics probed with the nuclear spin-lattice re- 
laxation rate (1/Ti). l/Ti of '^^P was measured at each 
resonance peak by the saturation-recovery method, and 
was uniquely determined by a single component in whole 
measured range. The inset of Fig. [3] shows T dependence 
of 1/TiT at low field ^qH ~ 0.6 T parallel and perpen- 
dicular to the c-axis. Below 1.5 K, 1/TiT as well as K 
along both directions becomes constant, indicative of the 
formation of a Fermi-liquid (FL) state of heavy electrons. 
In general, 1/Ti probes spin fluctuations perpendicular 
to applied H, and thus l/Ti in || c and H J. c are 



FIG. 3. (Color online) (Main panel)T dependence of low en- 
ergy spin fluctuations parallel and perpendicular to the c-axis 
at ~ 0.6 T evaluated with 1 /TiT measured in _ff ± c and -ff || c 
(See eq. [2]). The in-plane spin fluctuations are dominant at 
low T, suggesting that the spin dynamics also possess XY- 
type anisotropy. (Inset) T dependence of 1/T\T at 10.3 MHz 
(~ 0.6 T) for 7/ _L c and at 10.3 MHz (~ 0.6 T), 31.7 MHz 
(~ 1.8 T), and 107.2 MHz (~ 6.2 T) for H \\ c. {1/TiT)h\\c 
is independent of H up to 6.2 T. 

described as, 

{1/T^)h\\c = 2(/io7„y5^ |i/±(<?,C^rc.)P 

(X 2^2^ I ^^(q,^ - 0)p, and 

9 

il/Ti)H±c = (M07n)'^ [|ifc((Z,Wrcs)|' + (g, W,cs) P] 

oc [|5||(g,w ^ 0)p + |5x(g,c. ^ 0)p] . 

9 

(2) 

Here \X{uj)\ denotes the power spectral density of a 
time-dependent random variable X(t), and A is assumed 
to be q-independent due to the metallic state. From 
these equations, we can decompose spin fluctuations 
along each direction as shown in the main panel of Fig. 
El So 1*^-1- (^'"^ 0)p is dominant at low T, since 
(l/TiT)j/||j, is almost twice larger than {1/TiT)h±c- 
This indicates that the spin dynamics also possess XY- 
type anisotropy. The XK-typc spin fluctuations have the 
predominance of ferromagnetic (FM) correlations as in- 
ferred from the Korringa relation between (l/TiT)//||c 
and K± in low-T FL state, which is consistent with the 
previous '^^P-NMR result [5| and with the experimental 
facts that CeFePO is close to FM instability [3, Q. 

The evolution of the spin dynamics against H was 
investigated for both directions. Figure |4] shows T de- 
pendence of {1/TiT)h±c below 4 T (approximately the 
metamagnetic field, Hm) (a) and above 4 T (b). Al- 
though {1/TiT)h\\c does not depend on H up to 6.2 T 
as shown in the inset of FiglJl {1/TiT)h±c changes sig- 
nificantly by H as shown in FigsH] (a) and (b). H de- 
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FIG. 4. (Color online)r dependence of il/TiT)H±c below 4 T 
(a), and above 4 T (b). (c): H dependence of {1/TiT)h±c 
at low T. {1/TiT)h±c shows a distinct maximum at around 
Hm, indicating the enhancement of N{Ef) related to meta- 
magnetic anomaly, (d) H - T phase diagram defined by Tmax 
where {1/TiT)h±c shows a maximum. Non FL behavior char- 
acterized by continuous increase in {1/TiT)h±c with decreas- 
ing T (broken lines are shown in (b)) was observed in a narrow 
field region intervening between low-field paramagnetic (PM) 
state and high-field polarized PM state above 6 T. 

pendence of {1/TiT)h±c at 0.1 and 1.5 K is shown in 
Fig|4](c). {1/TiT)h±c shows a distinct maximum at Hm, 
suggesting that the enhancement of the density of states 
(DOS) is related to metamagnetic behavior. If we assume 
that 1/TiT cx NiEpf, N{Ef) at Hm is almost 1.5 times 
larger than N{Ey) at T. However, it is noteworthy that 
non FL behavior characterized by continuous increase in 
{1/TiT)ha_c with decreasing T was observed down to 

100 mK ((l/Tir)//j_c logT) in 4.3 T < ^ioH < 

5 T intervening between low-field paramagnetic (PM) 
state and high-field polarized PM state above 6 T. In the 
high-field polarized state, {1/TiT)h±c decreases with in- 
creasing H and {1/TiT)h±c at 14 T shows almost the 
same value as (LaCa)FcPO [~ 1.5 (sK)^^] without 4/ 
electrons 0. The H variation of (1/TiT)ha_c, reflecting 
the evolution of DOS at the Fermi level with H, strongly 
suggests the evolution of the Fermi surfaces (FSs) by H. 
It is noted that such a significant H dependence was not 
reported in the previous specific- heat measurement [sj . 
This would be because the magnetic field is applied in 
the various angles against the c-axis, and suggests that 
the metamagnetic behavior would be observed when H 
is exactly perpendicular to the c-axis. 

To investigate such evolution of FSs by H, we per- 
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FIG. 5. (Color online)Fermi surfaces (FSs) calculated with 
the ab initio band structure calculation as Ce-4/ electrons 
are itinerant (a), or localized (b), where the Fermi velocity 
is mapped on the FSs. (c) shows the FS colored by the jz 
character in the j — 5/2 multiplet of 4/ orbitals, where j 
is the total angular momentum, (d): the partial density of 
states. The Fermi level corresponds to eV. 

formed the ab initio band-structure calculation in the 
paramagnetic state of CeFePO by using the WIEN2k 
packagejlOj. FSs in the low-field region are composed 
of itinerant Ce 4/ electrons as shown in FiglS] (a). The 
large FS shows the characteristic neck structures around 
X-R at the Brillouin zone boundary, at which bound- 
ary electrons have small Fermi velocity or heavy electron 
mass. The orbital character is dominated by the jz — ± 
1/2 component in the j = 5/2 multiplet of 4/ orbitals, as 
seen in FigslS] (c) and (d) [ll|. These features of Fermi 
surface imply that the low-field HF state possesses the 
small q magnetic correlations and their in-plane compo- 
nent is much larger than the out-of- plane component, 
in good agreement with the experimental results. The 
band calculation also shows that applied field pinches off 
the neck FSs around R and X in order, that is, "f/ie 
field-induced Lifshitz transition^^ appears. This is accom- 
panied by a drastic change in DOS at the Fermi level, 
which can be a driving force for the metama gnet ic tran- 

With 
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sition and non-FL behavior around Hm 
applying higher H, 4/ electrons are localized, and thus 
the FSs become small. The resultant FSs are shown in 
FiglS] (b), familiar in the iron-based superconductors [l4| . 
Thus the scenario of the field-induced Lifshitz transition 
can link with the nature in the metamagnetic transition 
in CeFePO. 

Here, we compare the present results with CeRu2Si2, 
one of the most well-known metamagnetic compounds. 
Although both compounds show similar T and H de- 
pendence of Knight shift and l/TiT in H parallel to 
the magnetic easy axisfisj. as well as the similar H - 
T phase diagram defined by Tmax [13 as shown in Fig|3] 
(d), magnetic properties are quite different. For exam- 
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FIG. 6. (Color online) (a) Values of metamagnetic fields are 
plotted against Tmax determined from a maximum in static 
susceptibility (a) or inverse of Sommerfeld coefficient 7 at _ff 
= (b) in CeFePO, CeCue and CeRu2Si2 with pressurized 
and La-doped system. Broken lines are guide to eyes. A 
linear relation holds between two quantities, although three 
compounds possess quite different crystal structures and mag- 
netic properties, indicating that Hm is linked with the local 
Kondo singlet energy Tmax. 



pie, magnetic easy axis is different between CeFePO and 
CeRu2Si2: CeFePO possesses two dimensional Xy-type 
spin anisotropy, whereas CeRu2Si2 possesses Ising-type 
spin anisotropy [17|. As a result, the ground states of the 
crystal-field level are different and their metamagnetic 
behavior is observed in different directions. In addition, 
dominant magnetic fluctuations in CeFePO differ from 
those in CcRu2Si2. It is reported that CeRu2Si2 is lo- 
cated close to antiferromagnetic instability accompanied 
with FM fluctuations [l8l[ia|. 

Figure [S] shows the relationship between the metam- 
agnetic field Hm and the temperature where the bulk 
susceptibility shows a maximum Tmax or inverse of Som- 
merfeld coefficient 7 at = for CeCue, CeFePO. and 
CeRu2Si2 with doped and pressurized systems [sl. UtI. [20I- 
[2^. It deserves to note that the linear relation holds 
between the two quantities, notwithstanding that the 
three compounds possess totally different crystal struc- 
tures and magnetic properties. Since Tmax is regarded 
as a Kondo temperature Tk and approximately, the re- 
lation of jTk = const, holds in the HF state, these facts 
indicate that Hm is merely related to the local Kondo 
singlet energy Tmax and is not linked with the magnetic 
fluctuations originating from the intersite coupling be- 
tween neighboring Ce ions and/or the nesting between 
the "large" FS. The experimental facts that Hm is lin- 
early proportional to Tk in Figl6] strongly suggests that 
the metamagnetic behavior is linked with the Kondo 
breakdown 2^. Therefore, in the Ce-based metamagnets, 
the Kondo breakdown and the Fermi-surface instability 



accompanied by the drastic change of DOS occur almost 
simultaneously around Hm, which can induce novel non- 
FL behavior. 

In summary, wc performed "^^P-NMR in the uniaxially- 
aligned CcFcPO and found that CcFcPO possesses two 
dimensional XF-typc FM fluctuations, and shows meta- 
magnetic behavior when H is applied to H ± c below 
5 K, accompanied with non FL behavior around meta- 
magnetic field Hm ~ 4 T. As far as we know, this is a 
first example that the metamagnetic behavior occurs in a 
nonmagnetic Ce-based HF compound with the XY-type 
spin anisotropy. From the band calculation and the com- 
parison with other Ce-based metamagnets, we claim that 
Hm is a magnetic field breaking the local Kondo singlet, 
which is determined with the intrasite coupling between 
Ce-4/ and conduction electrons, and that the FSs change 
drastically due to the Kondo breakdown. 
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